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Abstract 



The Y(1S), Y(2S), and Y(3S) production cross sections are measured using a data 
sample corresponding to an integrated luminosity of 35.8 ± 1.4pb^^ of proton- 
proton collisions at ^/s = 7 TeV, collected with the CMS detector at the LHC. The 
Yresonances are identified through their decays to dimuons. Integrated over the Y 
transverse momentum range pj < 50GeV/c and rapidity range \y^\ < 2.4, and as- 
suming unpolarized Y production, the products of the Y production cross sections 
and dimuon branching fractions are 

£7(pp ^ Y(is)x) ■ i3(Y(is) ^ ^+^-) 

£7(pp ^ Y(2S)X) ■ i3(Y(2S) ^ ^+^-) 
a{pp Y(3S)X) ■ i3(Y(3S) ^ pl+^-) 

where the first uncertainty is statistical, the second is systematic, and the third is from 
the uncertainty in the integrated luminosity. The differential cross sections in bins of 
transverse momentum and rapidity, and the cross section ratios are presented. Cross 
section measurements performed within a restricted muon kinematic range and not 
corrected for acceptance are also provided. These latter measurements are indepen- 
dent of Y polarization assumptions. The results are compared to theoretical predic- 
tions and previous measurements. 



= (8.55 ± 0.05 toit ± 0-34) nb, 
= (2.21 ± 0.03 toil ± 0.09) nb, 
= (1.11 ± 0.02 toll ± 0-04) nb. 
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1 Introduction 

No existing theoretical approach successfully reproduces both the differential cross section and 
the polarization measurements of the ]/iporY states (l) in proton-proton collisions. Studying 
quarkonium hadroproduction at high center-of-mass energies and over a wide rapidity and 
transverse momentum range will facilitate significant improvements in our understanding of 
the processes involved. 

The first measurement of Y(nS) production, where n = 1-3, at\/s = 7TeV at the Large Hadron 
Collider (LHC) was reported by the Compact Muon Solenoid (CMS) Collaboration m, using 
a data sample corresponding to an integrated luminosity of 3pb^^. This paper constitutes an 
extension of that first cross section measurement, using a larger, independent sample, corre- 
sponding to an integrated luminosity of 35.8 ± 1.4 pb^^ collected in 2010. 

Two different approaches to the measurement of the Y production cross section are pursued 
in this paper. In each approach, the Y is reconstructed in the decay Y — > ■ In the first 

approach, a cross section measurement corrected for detector acceptance and efficiencies is 
presented, as in Ref. [2J. This cross section measurement depends on the spin alignment of the 
Y. No net polarization is assumed for the main results. To show the sensitivity of the results 
to the polarization and to allow for interpolation, we provide measurements for other polar- 
ization assumptions. Recently, the CMS Collaboration has measured the polarizations of the 
Y(nS) in pp collisions at ^/s = 7TeV, which are found to be small IS]. Cross section measure- 
ments are also provided in the Y transverse momentum (pj) and rapidity (y^) ranges matching 
those of the polarization measurement, and these polarization results are used to estimate the 
associated systematic uncertainty. The motivation for the second approach, also used by the 
ATLAS Collaboration f?|, is to eliminate the dependence of the measured cross sections on the 
spin alignment of the Y. In this second approach, a fiducial cross section measurement, cor- 
rected for detector efficiencies but not for acceptance, is presented. This cross section is defined 
within a muon kinematic range. 

The paper is organized as follows. Section |2] contains a short description of the CMS detector. 
Section |3] presents the data collection, the trigger and offline event selections, and the recon- 
struction of the Y resonances. Section |4] describes the measurement technique. The detector 
acceptance and efficiencies to reconstruct Y resonances that decay to two muons are discussed 
in Sections [Sjand [6} The evaluation of systematic uncertainties in the measurements is described 
in Section [yj^ In Sections [s] and |9j the Y(nS) fiducial and acceptance-corrected cross section re- 
sults and comparisons to other experiments and to theoretical predictions are presented. 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m inner diam- 
eter, producing a magnetic field of 3.8 T. Within the superconducting solenoid volume are 
a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a 
brass/ scintillator hadron calorimeter. Muons are detected by three types of gas-ionization de- 
tectors embedded in the magnet steel return yoke surrounding the solenoid: drift tubes, cath- 
ode strip chambers, and resistive-plate chambers. The muon measurement covers the pseudo- 
rapidity range |?/''| < 2.4, where = — ln[tan(0/2)] and the polar angle 9 is measured from 
the axis pointing along the counterclockwise-beam direction. The muon transverse momen- 
tum measurement, pj, based on information from the silicon tracker alone, has a resolution of 
about 1% for a typical muon in this analysis. The two-level CMS trigger system selects events 
of interest for permanent storage. The first trigger level, composed of custom hardware pro- 
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3 Data selection and event reconstruction 



cessors, uses information from the calorimeter and muon detectors to select events in less than 
3.2 }is. The high-level trigger software algorithms, executed on a farm of commercial proces- 
sors, further reduce the event rate using information from all detector subsystems. A detailed 
description of the CMS detector can be found in Ref. |,5J. 



Data are included in the analysis for all periods where the silicon tracker, the muon detec- 
tors, and the trigger were performing well and the luminosity information was available. The 
data were collected with two triggers. The trigger used in the first data-taking period requires 
the detection of two muons without an explicit pj requirement. The minimum distance be- 
tween each reconstructed muon trajectory and the average proton-proton interaction point 
in the transverse plane must be less than 2 cm. As the LHC instantaneous luminosity in- 
creased, this trigger was prescaled; consequently, a second trigger was used instead, which 
requires the two muons to have opposite charge and an invariant mass in the mass range 
1.5 < Mpifi < 14.5GeV/c^, in addition to the requirements imposed by the first trigger. All 
three muon systems take part in the trigger decision. The trigger used in the first data-taking 
period selected about 2 million events and the second trigger selected about 5 million events. 

Simulation is employed to design the offline selection, assess the detector acceptance, and study 
systematic effects. The Y(nS) events are simulated using pythia 6.412 |6|, which generates 
events based on the leading-order color-singlet and color-octet mechanisms, with nonrelativis- 
tic quantum chromodynamics (QCD) matrix elements, tuned by comparing calculations with 
CDF data [7], and applying the normalization and wave functions recommended in Ref. HI. 
The underlying-event simulation uses the CTEQ6L1 parton distribution functions [9]. Since 
PYTHIA does not provide a simulation of Y(2S) and Y(3S), the predictions for these states 
are obtained by replacing the Y(1S) mass in the simulation with the Y(2S) and Y(3S) masses, 
respectively. For the Y(2S) feed-down component, the masses of the 2P states replace the corre- 
sponding IP states. For the Y(3S) the feed-down is expected to be small and is not simulated. 
Final-state radiation (FSR) is implemented using PHOTOS [10. 11]. The response of the CMS 
detector is simulated with a GEANT4-based [12J Monte Carlo (MC) simulation program. Sim- 
ulated events are processed with the same reconstruction algorithms used for data. 

The offline selection starts from Y candidates reconstructed from pairs of oppositely charged 
muons with invariant mass between 7 and 14GeV/c^. The muons are required to have one or 
more reconstructed track segments in the muon systems that are well matched to the extrap- 
olated position of a track reconstructed in the silicon tracker. Quality criteria are applied to 
the tracks to reject muons from kaon and pion decays. Tracks are required to have at least 11 
hits in the silicon tracker, at least one of which must be in the pixel detector, and a track-fit 
per degree of freedom smaller than 5. In addition, tracks are required to extrapolate back to a 
cylindrical volume of radius 2 mm and length 25 cm, centered on the pp interaction region and 
parallel to the beam line. After offline confirmation of the trigger selection, muons are required 
to satisfy a kinematic threshold that depends on pseudorapidity 



3 Data selection and event reconstruction 



p'j > 3.75GeV/c 

> 3.5GeV/c 
pl^ > 3.0GeV/c 



if \r]i^'\< 0.8, 

if 0.8 < < 1.6, 

if 1.6 < < 2.4. 



(1) 



These kinematic acceptance criteria are chosen to ensure that the trigger and muon reconstruc- 
tion efficiencies are high and not rapidly changing within the phase space of the analysis. The 
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longitudinal separation between the two muons along the beam axis is required to be less than 
2 cm. The two muon helices are fit with a common vertex constraint, and events are retained if 
the fit "x^ probability is larger than 0.1%. If multiple dimuon candidates are found in the same 
event, the candidate with the smallest vertex-fit probability is retained; from simulation, the 
fraction of Y decays rejected by this requirement is about 0.1%. After offline selection, the data 
sample consists of approximately 573 000 events. 



4 Measurement of the inclusive differential cross section 

The product of the Y(nS) differential cross section, u, and the dimuon branching fraction, B, 
is determined from the signal yield ^'^(^^y corrected by the acceptance A and the efficiency e, 
using 

dc.(pp^Y(nS)X) _ _ W^^,p\,f;Ae) 

di/Y ^ ^ ^ ^ ^' F j - £ . . ' 

where L is the integrated luminosity of the data set, and Apj and Ay^ are the bin widths of the 

Y transverse momentum and rapidity, respectively. The rapidity is defined as y = j ln( g^|^^^ ), 
where £ is the energy and is the momentum component parallel to the beam axis of the 
muon pair. 

The Y(nS) yields are extracted via an extended unbinned maximum-likelihood fit to the dimuon 
invariant-mass spectrum. The measured mass line shape of each Y state is parametrized by a 
"Crystal Ball" (CB) [13J function, which consists of a Gaussian core portion and a power-law 
low-side tail to allow for FSR, with the low-mass tail parameters fixed from MC simulation [2J. 
The three Y(nS) states are fitted simultaneously since the three resonances overlap in the mea- 
sured dimuon mass range. The resolution, given by the standard deviation of the Gaussian 
component of the CB, is a free parameter in the fit, but is constrained to scale with the ratios 
of the resonance masses. However, the mass resolution varies with Y rapidity. Consequently, 
a single resolution term in the Gaussian component of the CB is not sufficient to describe the 
data. For this reason, in the pj intervals with sufficient statistical precision, the sum of two 
CBs with the same mean and FSR tail parameters, but different resolutions, is used for each 

Y state. The Y(nS) mass ratios are fixed to their world-average values [14J. The background in 
the 7-14 GeV/ (p- mass-fit range is nonpeaking and in some kinematic bins has a turn-on caused 
by the trigger and offline requirements. In general, the product of an error function and an 
exponential is chosen to describe the background [15 1, except when, for bins with poor statis- 
tical precision, a single exponential function is used. The dimuon invariant-mass spectra in 
the Y(nS) region, before accounting for acceptance and efficiencies, are shown in Fig. [ijand in 
Appendix [A| 

Following Ref. fZ^, given the significant Pj and dependencies of the acceptances and efficien- 
cies of the muons from Y(nS) decays, we correct for them on a candidate-by-candidate basis 
before performing the mass fit to obtain ^^IvS) ^^^^^ i^i Eq. (j2|. The fiducial differential cross 
section is determined from the efficiency-corrected signal yield within the kinematic region 
defined in Eq. ([ij. 



5 Acceptance 

The Y — > V-^Y-^ acceptance of the CMS detector is the product of two terms. The first is, for 
a given pj and y^, the fraction of dimuon decays in which both muons are within the phase 
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5 Acceptance 
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Figure 1: The dimuon invariant-mass distribution in the vicinity of the Y(nS) resonances for 
\y^\ < 2.4 (left) and for the subset of events where the rapidity of the Y(nS) satisfies \y^\ < 0.4 
(right). The solid lines represent the results of the fits to the signal-plus-background functions 
described in the text. 

space specified in Eq. The second is the probability that when there are only two muons in 
the event both can be reconstructed in the tracker without requiring the quality criteria. Both 
components are evaluated by simulation and parametrized as a function of pj and . 

Following Ref. |2J, the acceptance is defined by the ratio 



Pj, y^ I Si tracks satisfying Eq. ([T 



(3) 



where p'j, y'^ and Pj, y^ are the generated and reconstructed values, respectively, thus ac- 
counting for the detector resolution. The parameter N^'^" is the number of generated Y in a 
(p'j,i/'^) bin, while is the number of Y reconstructed with silicon tracks satisfying Eq. ([ij 
in the corresponding (pj, y^) bin. In addition the numerator requires the two tracks to be re- 
constructed with opposite charges and have an invariant mass within the Y mass-fit range of 
7-14GeV/c2. 

The acceptance is evaluated with a signal MC simulation sample in which the Y decay to two 
muons is generated with the EVTGEN |16| package, including FSR. There are no particles in 
the event besides the Y, its daughter muons, and the FSR photons. The Y mesons are gener- 
ated uniformly in pj and y^ . This sample is then simulated and reconstructed with the CMS 
detector simulation software to assess the effects of multiple scattering and finite resolution of 
the detector. An acceptance map with the assumption of zero Y polarization can be found in 
Ref. [2J. Systematic uncertainties arising from the dependence of the cross section measure- 
ment on the MC simulation description of the pj spectrum and resolution are evaluated in 
Section|7| The acceptance is calculated as a two-dimensional grid in pj and \y^\ using bin sizes 
of 0.1 in rapidity and 0.5 GeV/c in p^ for < p^ < 2 GeV/c and 1 GeV/c for 2 < p^ < 50 GeV/c. 
The corresponding correction is then performed on a candidate-by-candidate basis. The accep- 
tance varies with the Y(nS) mass. It decreases with rapidity, and there are no accepted events 
beyond | = 2.4. The acceptance has a minimum near pj = 5 GeV/c, as a result of the softer 
muon failing the Pj cut. The polarization of the Y strongly influences the muon angular dis- 
tributions and could be a function of pj. In order to show the sensitivity of the result to the 
Y(nS) polarization and to allow for interpolation, we provide cross section measurements for 
unpolarized (default) and 6 polarization scenarios in which the polar anisotropy parameter Ag 
is changed from fully longitudinal to fully transverse polarization, corresponding to = — 1, 
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—0.5, —0.25, 0.25, 0.5, 1, in both the center-of-mass helicity and ColHns-Soper flTl reference 
frames. Cross section measurements for the pj and ranges used in Ref. [3] are also pro- 
vided. In that case, the polarization results from Ref. [3] are used to estimate the corresponding 
systematic uncertainty. 



6 Efficiency 

The total muon efficiency is factorized into the three conditional terms, 

e = £(trig|id) x e(id|track) x e (track | accepted) = £trig x gj^ x Cfrack- (4) 

The tracking efficiency, £track/ combines the efficiency that the accepted track of a muon from 
a Y(nS) decay is reconstructed in the presence of additional particles in the silicon tracker, 
as determined with a track-embedding technique [18| , and the efficiency for the track to sat- 
isfy the track-quality criteria. The efficiency of the track-quality criteria, measured using a 
tag-and-probe (T&P) method on data IllSl , is nearly uniform and has an average value of 
(98.66 ± 0.05)%. A ]/xp data sample is utilized for the T&P efficiency measurement as it pro- 
vides a statistically independent, large-yield dimuon sample. The implementation of T&P for 
this analysis follows Ref. m. The muon identification efficiency, Eid, is the probability that the 
silicon track caused by a muon is correctly identified as a muon. The efficiency that an iden- 
tified muon satisfies the trigger is denoted by etrig. The muon trigger and muon identification 
efficiencies are also determined using the T&P technique. 

The Y efficiency is estimated from the product of the single-muon efficiencies. A factor, p, is 
used as a correction to this factorization hypothesis, and to account for possible biases intro- 
duced by the T&P efficiency measurement with the ]/xp sample. We define p as 

where e(Y) is the efficiency for a Y to pass the trigger and muon identification selections, and 
^(i" j/i/;) the corresponding efficiencies for positively and negatively charged 
muons from a J/t/? decay with the same pj and ?/ as a muon in the Y decay. The Y efficiency 
is taken from MC simulation generator-level matching, which is performed by associating the 
two generated muons from the Y with the reconstructed muons or trigger objects. The single- 
muon efficiencies are from the T&P method utilizing a]/ip MC simulation sample. Finally, the 
efficiency of the vertex-fit probability requirement is determined from data to be (99.16 ± 
0.09)% and constant over the entire kinematic range. 



7 Systematic uncertainties 

Systematic uncertainties in the cross section measurement stem from variations in the accep- 
tance determination, potential residual inaccuracies in the efficiency measurement, the method 
of yield extraction, and the integrated luminosity. For each uncertainty, we give below in paren- 
theses a representative range of values corresponding to the variation with pj. The acceptance 
is varied in the dimuon invariant-mass fit coherently by ±1 standard deviation, reflecting the 
uncertainty from the finite MC simulation statistics (0.3-1%). The acceptance is sensitive to 
biases in track momentum and differences in resolution between simulation and data. To de- 
termine the effect on the Y acceptance, we introduce a track pj bias of 0.2%, chosen based on 
the momentum scale biases seen in simulation and data [19]. We also vary the transverse mo- 
mentum resolution by ±10%, corresponding to the uncertainty in the resolution measurement 
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7 Systematic uncertainties 



using ]/xp in data. This reflects a conservative estimation of resolution effects. The acceptance 
map as a function of pj and \y^\ is then recalculated, and the systematic uncertainty is the 
difference in the resulting cross sections when using the perturbed acceptance map rather than 
the nominal one (0.0-0.7%). Imperfect knowledge of the production pj spectrum of the Y reso- 
nances at a/s = 7 TeV contributes a systematic uncertainty. Using either a flat pj distribution or 
the Pj distribution from PYTHIA, which is found to be consistent with the previously measured 
Pj distribution [2J, gives rise to a systematic uncertainty (0.2%). FSR is incorporated into the 
simulation using the PHOTOS algorithm. To estimate the systematic uncertainty associated with 
this procedure, the acceptance is calculated without FSR, and 20% of the difference is taken as 
the uncertainty (0.1-0.8%), based on a study in Ref. |[TT|. 

Variation of the measured factorized efficiencies within their uncertainties also gives rise to 
a systematic uncertainty. The systematic uncertainties for the tracking efficiency (0.3-0.4%), 
muon identification efficiency (2-4%), and trigger efficiency (1-5%) are evaluated conserva- 
tively by coherently varying all bins by ±1 standard deviation. The systematic uncertainty 
arising from the choice of bin size for the efficiencies is determined by fitting the efficiency 
turn-on curves as a function of muon pj in different regions using a hyperbolic tangent 
function and taking the muon efficiencies from the function instead of the binned value to 
compute the cross section (1-4%). The intrinsic bias from the T&P method, including possible 
bias in the T&P technique and differences in the ]/xp and Y kinematics, as well as the possible 
misestimation of the double-muon Y efficiency as the product of the single-muon efficiencies, 
are all included in the correction factor p. As a conservative estimate of the systematic uncer- 
tainty associated with p, the measurements are repeated with a correction factor of unity and 
half of the variation is taken as the systematic uncertainty (2-5%). 

In addition, systematic uncertainties may arise from differences between the dimuon invariant- 
mass distribution in the data and the probability density functions (PDF) chosen for the signal 
and background components in the fit. Since the CB parameters, which describe the radiative 
tail of each signal resonance, are fixed from MC simulation in the fit to the data, we fit the 
full data set with free tail parameters and use the values obtained to fix the tail parameters 
for the yield extraction in the (Apj, Ay^) bins. The difference in the fit yield is taken as a 
systematic uncertainty (1-4%). We vary the background PDF by replacing the exponential 
convolved with an error function by a polynomial function, while restricting the fit to the mass 
range 8-12GeV/c^ (1-5%). The determination of the integrated luminosity is made with an 
uncertainty of 4% [20]. A summary of systematic uncertainties for the Y(nS) production cross 
section, integrated over the full transverse momentum (pj) and rapidity (y^) ranges, is shown 
in Table [l] The largest sources of systematic uncertainty arise from the statistical precision of 
the efficiency measurements determined from data, the efficiency correction factor p, and from 
the measurement of the integrated luminosity. 

The cross section measurement uses acceptance maps corresponding to different Y polarization 
scenarios. The values of the resulting cross sections vary approximately linearly by about ±5%, 
±10%, and ±20%, respectively, assuming \g = ±0.25, ±0.5, and ±1, as shown in Table|2j The 
cross sections are also measured for 10 < pj < 50GeV/c and \y^\ < 1.2 in which we use the 
measured Y(nS) polarizations [3] to compute the acceptance corrections. The three anisotropy 
parameters in the center-of-mass helicity and Collins-Soper frames are varied coherently by 
±1 standard deviation, and the largest positive and negative variations with respect to the 
nominal (no polarization) case are taken as systematic uncertainties. These are listed in Table |4] 
They are comparable to, or smaller than, the result of varying the longitudinal or transverse 
polarizations by setting Aq to ±0.25 for the Y(1S) case, while they are between the results 
obtained by setting Aq to ±0.25 and ±0.5 for the Y(2S) and Y(3S). The fiducial cross sections 
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do not depend on the acceptance, the assumed Y polatization, or the associated uncertainties. 
The definition of the acceptance in Eq. (|3| includes reconstructed quantities. The variation in 
the cross section using only generator-level quantities is less than 1%. 

Table 1: Relative systematic uncertainties in the Y(nS) production cross section, integrated over 
the rapidity range \y^\ < 2.4, times the dimuon branching fraction, in percent. The symbols 
A, ex&p, £p, and PDF refer to the systematic uncertainties arising from the acceptance, tag-and- 
probe efficiencies, correction factor p, and signal-and-background PDF. The remaining system- 
atic uncertainties are summed in the "other" category. The integrated luminosity uncertainty 
of 4% is not shown. The numbers in parentheses are negative variations. 





Pt (GeV/c) 


A 


^T&P 


£p 


PDF 


other 


Y(1S) 


0-50 


1.0(1.0) 


5.2 (4.3) 


3.4 


1.8 


0.4 (0.3) 


Y(2S) 


0^2 


1.1 (1.1) 


5.5 (4.1) 


3.7 


2.6 


0.4 (0.4) 


Y(3S) 


0-38 


1.2(1.1) 


6.7(4.9) 


4.0 


3.8 


0.6 (0.5) 



Table 2: The fractional change in percent to the central value of the Y(nS) cross section inte- 
grated over the rapidity range \y^\ < 2.4, relative to the unpolarized value, for six polarization 
scenarios in the center-of-mass helicity and Collins-Soper frames. The polarization assumption 
changes from fully longitudinal to fully transverse polarization as Ag changes from —1 to 1. 





Helicity Frame \g 


Collins-Soper Frame \g 




Pt (GeV/c) 


1 


0.5 


0.25 -0.25 


-0.5 


-1 


1 


0.5 0.25 -0.25 -0.5 


-1 


Y(1S) 


0-50 


+19 


+10 


+5 -5 


-11 


-24 


+16 


+8 +4 -5 -9 


-19 


Y(2S) 


0-42 


+14 


+5 


+3 -7 


-12 


-24 


+13 


+6 +2 -6 -10 


-20 


Y(3S) 


0-38 


+16 


+9 


+5 -4 


-9 


-21 


+14 


+8 +5 -3 -7 


-17 



8 Differential fiducial cross section measurement and compari- 
son to theory 

The fiducial Y(nS) cross sections are determined from the efficiency-corrected signal yields 
within the muon kinematic range specified by Eq. ([l|, using Eq. ^ with the acceptance term 
set to unity. The resulting total fiducial Y(nS) cross sections times dimuon branching fractions 
for \y^\ < 2.4 are 

(7(pp ^ Y(IS)X) ■ 6(Y(1S) ^+|^-) = (3.06 ± 0.02 tofg ± 0.12) nb, 
(7(pp ^ Y(2S)X) ■ 6(Y(2S) -^ll+^-) = (0.910 ± 0.011 ± 0.036) nb, 
(7(pp ^ Y(3S)X) ■ 6(Y(3S) ^ ^^^-) = (0.490 ± 0.010 ±0.020) nb, 

where the first uncertainty is statistical, the second is systematic, and the third is associated 
with the estimation of the integrated luminosity of the data sample. The integrated results 
are obtained from the sum of the differential pj results. The measured cross sections include 
feed-down from higher-mass bottomonium states. 

The Y(nS) differential px fiducial cross sections are plotted in Fig.[2]and summarized in Table|3] 
In the figures, i3(Y(nS) — > is denoted as B(/^f/). The results are given for six rapidity 

intervals in Fig. |3] Here, and throughout this paper, in figures illustrating differential cross 
sections, the data points are plotted at the average of the px (or rapidity) of the data in each 
bin. The pj dependence of the cross sections has the same trend for all six rapidity intervals. 
The Y(nS) px-integrated, differential rapidity fiducial cross sections are plotted in Fig. Ill The 
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8 Differential fiducial cross section measurement and comparison to theory 



fiducial cross sections are all roughly constant from \y^\ = to about 1.6, where they then fall 
quickly. The ratios of the Y(nS) differential px fiducial cross sections are shown in Fig. [s] The 
ratios increase with pj. The Y(3S) is not expected to have a significant feed-down component 
from decays of higher bottomonium states, unlike the Y(2S) and Y(1S). Therefore, a possible 
explanation for the px dependence of the ratios is that the px dependence of the cross section 
for direct production of each Y resonance is the same, and the softening of the Y(1S) and Y(2S) 
px spectra is due to feed-down from higher-mass states. 
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Figure 2: Measured Y(nS) differential fiducial cross sections as a function of pj in the rapidity 
range \y^\ < 2.4. The error bars indicate the total uncertainties, except for the 4% uncertainty 
in the integrated luminosity. 

A comparison between the fiducial cross section measurement and theoretical predictions is 
shown in Fig. [6] Each of the predictions is for unpolarized Y(nS) production. The compar- 
ison is made to the CASCADE [21J MC generator in the fixed-order-plus-next-to-leading-log 
(FONLL) framework, including feed-down from ;\;b(lP), ;\;j,(2P), Xbi^^) [22J, and other higher- 
mass Y states, and to pythia [8] including feed-down for the Y(1S) and Y(2S) from the P-wave 
states with the same principal quantum number. The px dependence of the cross section pre- 
dicted by CASCADE agrees with the data for the Y(1S), is marginally consistent for the Y(2S) 
but does not describe the Y(3S) spectrum, where it predicts a softer px spectrum. The CAS- 
CADE predictions also describe the rapidity dependence of the Y(1S) measurement, but not of 
the Y(2S) or Y(3S) as shown in Fig.|6](d). For each resonance, the total cross section predicted 
by PYTHIA is higher, by factors of about 2, than the measured cross section. In Fig. |6j for each 
resonance the pythia prediction is normalized to the measured total cross section, in order to 
facilitate the comparison of the cross section px dependence with the prediction. The pythia 
prediction of the px dependence agrees with data for the Y (IS) and Y(3S),butnot for the Y(2S); 
PYTHIA provides a good description of the rapidity dependence of the three states. Complete 
tables of results for the differential fiducial cross sections for the three Y states are available in 
Appendix [A| 
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Figure 3: The Y(nS) differential fiducial cross sections as a function of pj in six rapidity ranges, 
scaled by the factors shown in the figures for clarity. 



10 9 Acceptance-corrected differential cross section measurement and comparison to theory 



J3 
C 

3 
pa 

X 
T3 



10' 



CMS, \fs = 7 TeV 
L = 36 pb-\ < 50 GeV/c 



o Y(1S) 
□ Y(2S) 
A Y(3S) 

Fiducial 

I I I I I I I 
0.5 1 



I Stat. unc. 
total unc. except 
luminosity (4%) 



I I I 
1.5 



|y^ 



Figure 4: Differential fiducial cross sections of the Y(nS) as a function of rapidity. The regions 
show the statistical uncertainties and the error bars show the total uncertainties, except for the 
4% uncertainty in the integrated luminosity. 
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9 Acceptance-corrected differential cross section measurement 
and comparison to theory 

The acceptance-corrected Y(nS) production cross sections times the dimuon branching frac- 
tions at ^/s = 7TeV for \y^\ < 2.4 and pj < 50GeV/c are measured to be 



C7(pp ^ Y(IS)X) ■ e(Y(lS) 
cr(pp ^ Y(2S)X) ■ S(Y(2S) 
cr(pp ^ Y(3S)X) ■ i3(Y(3S) 



(8.55 ± 0.05 Zq^ ± 0.34) nb, 
(2.21 ± 0.03 ± 0.09) nb, 
(1.11 ±0.02 ± 0.04) nb. 



(6) 



where the first uncertainty is statistical, the second is systematic, and the third is from the 
estimation of the integrated luminosity. These results assume unpolarized Y(nS) production. 
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Figure 6: Differential fiducial cross section of (a) Y(1S), (b) Y(2S), and (c) Y(3S) as a function 
of pj in the rapidity range \y^\ < 2.4, and comparison to the predictions from CASCADE and 
PYTHIA. (d) Differential fiducial cross section of the Y(nS) as a function of rapidity and compar- 
ison to the predictions from CASCADE and PYTHIA. The PYTHIA prediction is normalized to the 
measured total cross section. The bands indicate the estimated uncertainties in the CASCADE 
prediction. 
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Table 3: The product of the fiducial or acceptance-corrected Y(nS) production cross sections, a, 
integrated and differential in pj, and the respective dimuon branching fraction, B, integrated 
over the rapidity range \y^\ < 2.4. The cross sections assume the Y(nS) are unpolarized. The 
fiducial Y(nS) cross sections are independent of the Y(nS) polarization. The statistical un- 
certainty (stat.), the sum of the systematic uncertainties in quadrature (Egyst.)/ ^^id the total 
uncertainty (Aa; including stat., Dgyst./ ^rid the uncertainty in the integrated luminosity) are in 
percent. The numbers in parentheses are negative variations. 



Fiducial Cross Section 





Pt (GeV/c) 


mean 


a ■ B (nb) 


stat. 


1-isyst. 

IT 




a ■ B (nb) 


stat. 


i-isyst. 
IT 






0-0.5 


0.33 


0.0440 


5.4 


8(8) 


11(11) 


0.0859 


5.4 


8(7) 


11 (10) 




0.5-1 


0.77 


0.133 


3.1 


8(8) 


10(10) 


0.263 


3.3 


8(7) 


9(9) 




1-1.5 


1.26 


0.182 


2.5 


8(8) 


9(9) 


0.374 


2.6 


8(8) 


9(9) 




1.5-2 


1.75 


0.228 


2.4 


8(8) 


10 (9) 


0.505 


2.4 


9(8) 


10 (9) 




2-3 


2.49 


0.442 


1.6 


8(7) 


9(8) 


1.16 


1.6 


8(10) 


9(11) 




3-4 


3.48 
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1.8 


6(6) 


8(7) 


1.21 


2.1 


7(6) 


9(8) 
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4.48 
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1.8 


7(7) 


8(8) 
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Cross Section 



The Y(1S) integrated production cross section in the restricted rapidity range \y^\ < 2.0 is 
7.496 ± 0.052 (stat.) nb, which is consistent with the previous CMS result of 7.37 ± 0.13 (stat.) nb 
121, measured in the same rapidity range. The results of the Y(nS) production cross sections for 
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the same pj and ranges used for the measurement of the Y(nS) polarizations in Ref. ||3l are 
shown in Table H) 

Table 4: The product of the acceptance-corrected Y(nS) production cross sections, a, and the 
dimuon branching fraction, B, integrated over the rapidity range \y^\ < 1.2, and the pj range 
from 10 to 50GeV/c, as used in Ref. [3J for the measurement of the Y polarizations. The 
cross sections assume the Y(nS) are unpolarized. The statistical uncertainty (stat.), the sum 
in quadrature of the systematic uncertainties (Egyst.)/ excluding the contribution from the po- 
larization uncertainty, the systematic uncertainties from the polarization (pol.), and the total 
uncertainty (Aa; including stat., I^gygt./ pol., and the uncertainty in the integrated luminosity) 
are in percent. The numbers in parentheses are negative variations. 





a ■ B (nb) 


stat. 

(7 


^syst. 


pol. 


ha 


Y(1S) 
Y(2S) 
Y(3S) 


0.558 
0.213 
0.127 


1.3 
2.4 
3.2 


6(5) 
5(5) 
7(5) 


4(2) 
7(3) 
7(3) 


8(7) 

10 (8) 

11 (8) 



The acceptance-corrected Y(nS) differential px cross sections for the rapidity range \y^\ < 2.4 
are plotted in Fig. [7| and summarized in Table |3] Figure |8] shows the same for the ranges 
10 < pj < 50GeV/c, ly^l < 1.2 used in Ref. [3J and includes the systematic uncertainties from 
the polarization measurement of Ref. [3J, as explained in Section [t] The Y(nS) differential pj 
cross sections for six different rapidity bins are given in Fig.|9] The pj dependence of the cross 
section in the six exclusive rapidity intervals shows a similar trend within the uncertainties. The 



Y(nS) pT"iritegrated, differential rapidity cross section results are shown in Fig. 10 (a). Simi- 
lar to the fiducial differential rapidity cross sections, the acceptance-corrected cross sections 
are approximately flat from \y^\ = to about 2.0, where they then begin to fall. Comparison 



with similar results from the LHCb Collaboration |23| is shown in Fig. 10 (b). The two sets of 
measurements are complementary in their rapidity coverage and consistent within the uncer- 
tainties in the region of overlap. The fiducial cross sections and the acceptance-corrected cross 
sections exhibit similar pj and \y^\ dependencies. However, the decrease in the cross section 
at large values of the rapidity is greater for the fiducial cross section than for the acceptance- 
corrected cross section because the acceptance also decreases with rapidity. The ratios of the 
Y(nS) differential pj cross sections are plotted in Fig.jllJ along with comparisons to the CAS- 
CADE and PYTHIA predictions. The ratios increase with pj, as they did for the fiducial cross 
sections. The CASCADE prediction is consistent with the Y(2S)/Y(1S) and Y(3S)/Y(2S) mea- 
surements, while it disagrees with the Y(3S)/Y(1S) results at low pj. The pythia prediction 
agrees with the measured Y(3S)/Y(1S) values, but is inconsistent with the Y(2S)/Y(1S) and 
Y(3S)/Y(2S) results. 

The acceptance-corrected differential px and rapidity Y(nS) cross sections and the theoretical 



predictions are shown in Fig. 12 The measurements and predictions in Figs. 12 (a), 12 (b), and 
12 (c) are for \y^\ < 2.0 and assume unpolarized Y(nS) production. Comparisons are made 
to the CASCADE MC generator; the normalized PYTHIA (as explained in Section |8]l; the color- 
evaporation model (CEM) Il24l with feed-down not included; nonrelativistic QCD (NRQCD) to 
next-to-leading order (NLO) including feed-down, as described in Ref. f25|; the color-singlet 
model (CSM) to NLO and NNLO* [26|, with feed-down accounted for by scaling the Y(1S) 
and Y(2S) direct-production cross sections by factors 2 and 1.43, respectively [261 , and no 
feed-down for the Y(3S). The theoretical predictions are based on published models for un- 
polarized Y(nS) production, and, except for NRQCD |25l, are made for lower ^/s pTll24ll26l . 
These models have been updated by the authors to \/s = 7 TeV when relevant. The updates 
are unpublished and are in the form of private communications. Our measured Y(1S) cross 
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Figure 8: Acceptance-corrected differential cross sections as a function of pj for \y^\ < 1.2. 
The error bars represent the total uncertainties, including the systematic uncertainties from 
the measurement of the Y(nS) polarization [3J, but not the uncertainty (4%) in the integrated 
luminosity. 

section is smaller than the CEM prediction at low pj. The data agree with NRQCD above pj 
= lOGeV/c for the Y(1S) and with CASCADE for the Y(1S), but the agreement is not as satis- 
factory for the Y(2S) and Y(3S) at low momentum. The NLO CSM does not describe the data, 
while the NNLO* CSM shows better agreement within the large uncertainties. The total cross 
section predicted by PYTHIA is higher than the measured cross section, so in Fig.[l2]the PYTHIA 
prediction is normalized to the measured Y(nS) cross sections. The pj dependence of the cross 
section predicted by PYTHIA agrees with the data for the Y(1S) and Y(3S) but not for the Y(2S). 
PYTHIA also describes the rapidity dependence over the range of the measurement, as shown 



in Fig. 12 (d). Complete tables of results for the differential cross sections for the three Y states 
are available in Appendix |A| including variations for extreme polarization scenarios. 
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Figure 9: Differential production cross sections of the Y(nS) as a function of pj in six rapidity 
regions. The measurements have been multiplied by the factors shown in the figure for clarity. 

10 Summary 

Measurements of the Y(nS) differential and total production cross sections from proton-proton 
collisions at -^s = 7 TeV with the CMS detector have been presented. The results have been 
shown in two ways: as fully acceptance-corrected cross sections, and fiducial cross sections in 
which both muons from the Y(nS) decay are within the detector acceptance. The latter cross 
sections are independent of the assumed Y(nS) polarizations. The differential cross sections 
have been given as a function of pj and \y^\, and compared to theoretical predictions. The dif- 
ferential cross sections as a function of pj and for each Y(nS) state have also been measured 
and compared to theoretical predictions. Finally the Y cross section ratios have been given. 
The dominant sources of systematic uncertainty in the cross section measurements arise from 
the data-driven determination of the muon identification and trigger efficiencies, the correction 
factor for the efficiencies, and the integrated luminosity. 

The measurements are consistent with previous CMS results f2] based on less than 10% of the 
integrated luminosity analyzed here. These earlier measurements have been extended in terms 
of both the precision attained and the kinematic reach. In addition, this paper expands upon 
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Figure 10: (a) Acceptance-corrected differential production cross sections as a function of rapid- 
ity, and (b) comparison with LHCb results |23|. The bands represent the statistical uncertainty 
and the error bars represent the total uncertainty, except for those from the Y(nS) polarization 
and, in the case of (a), also the integrated luminosity. 

the previous result by the inclusion of fiducial cross section measurements and the polarization 
systematics, utilizing the recent Y polarization results from CMS |I3|. The results presented 
here will allow for a more precise determination of the parameters of the various bottomonium 
production models. 
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A Supplementary Material 

CMS,\fs = 7 TeV, L = 36 pb\ ly^ f < 2.4 





IL*IL invariant mass (GeV/c^) 

Figure A.l: Fit to the dimuon invariant-mass distribution in the specified px regions for \y\ < 
2.4, before accounting for acceptance and efficiency. The solid line shows the result of the fit 
described in the text, with the dashed line representing the background component. 



22 



A Supplementary Material 



CMS,\is = 7 TeV, L = 36 pis\ ly^ I < 2.4 




8 9 10 11 12 13 



8 9 10 11 12 13 



8 9 10 11 12 13 

invariant mass (GeV/c^) 



Figure A.2: Fit to the dimuon invariant-mass distribution in the specified px regions for \y\ < 
2.4, before accounting for acceptance and efficiency. The solid line shows the result of the fit 
described in the text, with the dashed line representing the background component. 
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Table A.l: The product of the fiducial Y(1S) production cross sections, a, integrated and 
differential in pj, and the dimuon branching fraction, B, integrated over the rapidity range 
ly"^! < 2.4. The statistical uncertainty (stat.), the sum of the systematic uncertainties in quadra- 
ture (Zsyst.)/ and the total uncertainty {Acr; including stat., Esyst./ and the imcertainty in the 
integrated Iiuninosity) are in percent. The numbers in parentheses are negative variations. 
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Table A.2: The product of the fiducial Y(2S) production cross sections, a, integrated and 
differential in pj, and the dimuon branching fraction, B, integrated over the rapidity range 
ly"^! < 2.4. The statistical uncertainty (stat.), the sum of the systematic uncertainties in quadra- 
ture (Zsyst.)/ and the total uncertainty {Acr; including stat., Esyst./ and the imcertainty in the 
integrated Iiuninosity) are in percent. The niunbers in parentheses are negative variations. 
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Table A.3: The product of the fiducial Y(3S) production cross sections, a, integrated and 
differential in pj, and the dimuon branching fraction, B, integrated over the rapidity range 
ly^l < 2.4. The statistical uncertainty (stat.), the siun of the systematic uncertainties in quadra- 
ture (Zsyst.)/ and the total uncertainty (A(r; including stat., Esyst' the imcertainty in the 
integrated limiinosity) are in percent. The numbers in parentheses are negative variations. 
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Table A.4: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in p^, and the respective dimuon branching fraction, B, integrated over the rapidity 



range [i/* | < 0.4. The statistical uncertainty (stat.), the sum of the systematic uncertainties in 
quadrature (Zsyst.)/ and the total uncertainty (Ac; including stat., I^gyst / ^ri'^ ^^e uncertainty in 
the integrated luminosity) are in percent. The numbers in parentheses are negative variations. 
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5.7 



12(11) 
11 (11) 
8(7) 
5(5) 
4(5) 



0-7 

7-12 
12-38 



0.0434 
0.0203 
0.0176 



6.5 
7.5 
5.9 



14 (12) 

10 (9) 
5(5) 



J2L. 



15 (14) 
11 (12) 
11 (11) 
15 (15) 
8(8) 
7(7) 
7(7) 



15 (14) 
12 (12) 
10 (10) 
9(9) 
8(9) 



16 (15) 
13 (13) 
9(9) 



Table A.5: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in pj, and the respective dimuon branching fraction, B, integrated over the rapidity 
range 0.4 < < 0.8. The statistical uncertainty (stat.), the sum of the systematic uncertain- 
ties in quadrature (Egygt.), and the total uncertainty (Acr; including stat., Egyst./ the uncer- 
tainty in the integrated luminosity) are in percent. The numbers in parentheses are negative 
variations. 

0.4 < ly^l < 0.8 



Y(1S) 



Y(2S) 



Y(3S) 



pT(GeV/c) 

0-2 

2-4 

4-6 

6-8 

8-11 

11-15 

15-50 



cr ■ B (nb) 
0.112 
0.127 
0.0967 
0.0695 
0.0731 
0.0467 
0.0314 



stat- 
er 
3.3 

3.0 
2.8 
3.3 
2.9 
3.3 
0.8 



lyst. 



12 fl2) 

9(9) 
8(8) 
12 (12) 
6(5) 
6(6) 
4(4) 



0-3 
3-7 
7-11 

11- 15 
15-42 

0-7 
7-12 

12- 38 



0.0472 
0.0580 
0.0334 
0.0154 
0.0147 
0.0504 
0.0247 
0.0155 



7.2 
4.5 
5.2 
6.6 
5.9 
6.8 
6.8 
6.5 



12(11) 
8(8) 
7(5) 
5(6) 
4(5) 

13(11) 
7(7) 
5(5) 



a 



13(13) 
10(11) 
9(9) 
13(13) 
8(7) 
8(8) 
6(6) 



14 (14) 
10(10) 

10 (8) 
9(10) 

8(8) 
15(13) 
11 (10) 

9(9) 
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Table A.6: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in pj, and the respective dimuon branching fraction, B, integrated over the rapidity 
range 0.8 < < 1.2. The statistical uncertainty (stat.), the sum of the systematic uncertain- 
ties in quadrature (Egyst.)/ and the total uncertainty [Act; including stat., I^gyst.' the uncer- 
tainty in the integrated liuninosity) are in percent. The niunbers in parentheses are negative 
variations. 



0.8 < ly^l < 1.2 





Pt (GeV/c) 


a ■ B (nb) 


stat. 




Act 




0-2 


0.107 


3.7 


9(9) 


11 (10) 




2-4 


0.160 


2.7 


9(9) 


10(10) 




4-6 


0.0995 


2.9 


7(6) 


8(8) 


Y(1S) 


6-8 


0.0770 


3.1 


12(12) 


13(13) 




8-11 


0.0746 


2.9 


5(8) 


7(9) 




11-15 


0.0464 


3.4 


6(6) 


8(8) 




15-50 


0.0313 


3.8 


5(6) 


7(8) 




0-3 


0.0561 


6.1 


11 (10) 


13(12) 




3-7 


0.0610 


5.0 


8(8) 


10(10) 


Y(2S) 


7-11 


0.0347 


5.3 


14 (14) 


15(15) 




11-15 


0.0168 


6.2 


6(7) 


9(10) 




15-42 


0.0142 


6.2 


6(7) 


10(10) 




0-7 


0.0526 


7.7 


17(17) 


19(19) 


Y(3S) 


7-12 


0.0264 


6.9 


11 (14) 


14(17) 




12-38 


0.0154 


6.8 


10(10) 


13 (13) 



Table A.7: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in pj, and the respective dimuon branching fraction, B, integrated over the rapidity 
range 1.2 < \y^\ < 1.6. The statistical uncertainty (stat.), the sum of the systematic uncertain- 
ties in quadrature (Zgyst.)/ and the total uncertainty {Acr; including stat., Esyst./ the uncer- 
tainty in the integrated liuninosity) are in percent. The niunbers in parentheses are negative 
variations. 



1.2 < ly^l < 1.6 





pT(GeV/c) 


(7 • B (nb) 


stat. 

(7 








0-2 


0.119 


3.5 


6(12) 


8(13) 




2-A 


0.184 


2.4 


7(6) 


8(8) 




4-6 


0.1059 


2.8 


6(5) 


8(7) 


Y(1S) 


6-8 


0.0747 


3.0 


5(6) 


7(7) 




8-11 


0.0697 


3.0 


5(5) 


7(7) 




11-15 


0.0456 


3.3 


6(6) 


8(8) 




15-50 


0.0319 


3.8 


7(8) 


9(10) 




0-3 


0.0582 


6.2 


6(6) 


10 (9) 




3-7 


0.0684 


4.5 


7(6) 


9(8) 


Y(2S) 


7-11 


0.0317 


5.5 


5(6) 


9(9) 




11-15 


0.0164 


6.5 


5(5) 


9(9) 




15^2 


0.0147 


6.4 


8(9) 


11 (12) 




0-7 


0.0589 


6.9 


6(5) 


10(10) 


Y(3S) 


7-12 


0.0214 


7.7 


9(9) 


13 (12) 




12-38 


0.0171 


6.7 


8(8) 


11(11) 



27 



Table A.8: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in pj, and the respective dimuon branching fraction, B, integrated over the rapidity 
range 1.6 < \y^\ < 2.0. The statistical uncertainty (stat.), the sum of the systematic uncertain- 
ties in quadrature (Egyst.), and the total uncertainty [Act; including stat., I^gyst.' the uncer- 
tainty in the integrated liuninosity) are in percent. The niunbers in parentheses are negative 
variations. 



1.6 < ly^l < 2.0 





Pt (GeV/c) 


cr ■ B (nb) 


stat. 

cr 




A(7 




0-2 


0.105 


3.1 


9(8) 


10 (9) 




2-4 


0.164 


2.5 


6(5) 


8(7) 




4-6 


0.1058 


2.7 


5(4) 


7(6) 


Y(1S) 


6-8 


0.0564 


3.2 


4(4) 


6(6) 




8-11 


0.0500 


3.4 


6(6) 


8(8) 




11-15 


0.0321 


4.1 


5(5) 


8(8) 




15-50 


0.0274 


4.4 


6(6) 


8(8) 




0-3 


0.0452 


6.2 


8(8) 


11 (11) 




3-7 


0.0661 


4.3 


7(7) 


9(9) 


Y(2S) 


7-11 


0.0205 


7.0 


10 (10) 


13(13) 




11-15 


0.0112 


8.3 


5(5) 


11 (11) 




15-42 


0.0120 


8.1 


6(5) 


11 (11) 




0-7 


0.0525 


7.2 


9(9) 


12(12) 


Y(3S) 


7-12 


0.0120 


10.6 


15 (15) 


19(19) 




12-38 


0.0124 


8.8 


6(6) 


11 (12) 



Table A.9: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in Pj, and the respective dimuon branching fraction, B, integrated over the rapidity 
range 2.0 < \y^\ < 2.4. The statistical uncertainty (stat.), the sum of the systematic uncertain- 
ties in quadrature (Egyst.)/ and the total uncertainty {Acr; including stat., Esyst./ and the uncer- 
tainty in the integrated liuninosity) are in percent. The niunbers in parentheses are negative 
variations. 





2.0 < \y^\ 


< 2.4 








pT(GeV/c) 


a ■ B (nb) 


stat. 


^syst. 






0-2 


0.0299 


5.3 


9(9) 


11 (11) 




2-4: 


0.0671 


2.8 


7(5) 


8(7) 




4-6 


0.0515 


3.7 


12 (12) 


13(13) 


Y(1S) 


6-8 


0.0216 


5.6 


6(7) 


9(10) 




8-11 


0.0194 


5.8 


6(6) 


9(9) 




11-15 


0.0113 


7.2 


9(8) 


12 (12) 




15-50 


0.0092 


9.1 


13 (12) 


16 (16) 




0-3 


0.0119 


13.2 


11 (10) 


18(17) 




3-7 


0.0278 


7.6 


9(8) 


12(12) 


Y(2S) 


7-11 


0.0098 


10.9 


9(8) 


14 (14) 




11-15 


0.00440 


14.7 


8(7) 


17 (17) 




15-42 


0.00365 


20.4 


10 (9) 


23 (23) 




0-7 


0.0223 


11.2 


10 (9) 


16(15) 


Y(3S) 


7-12 


0.0054 


19.3 


7(7) 


21 (21) 




12-38 


0.00518 


15.8 


8(8) 


18(18) 
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Table A.IO: The product of the fiducial Y(nS) production cross sections, a, integrated and dif- 
ferential in \y^\, and the respective dimuon branching fraction, B, integrated over the px range 
Pj < 50GeV/c. The statistical uncertainty (stat.), the sum of the systematic uncertainties in 
quadrature (Sgyst.), and the total uncertainty (Au; including stat., Yi,sysi.r the uncertainty in 
the integrated luminosity) are in percent. The numbers in parentheses are negative variations. 

< 50GeV/c 





1 1 

lyl 


a ■ B (nb) 


stat. 

(T 


a 


A(7 
(7 




0-0.2 


0.262 


1.7 


9 (9) 


10 (10) 




0.2-0.4 


0.268 


1.7 


9 (9) 


10 (10) 




0.4-0.6 


0.269 


1.7 


9 (9) 


10 (10) 


Y(1S) 


0.6-0.8 


0.295 


1.7 


7(6) 


8 (8) 




0.8-1.0 


0.295 


1.7 


8(8) 


9 (9) 




1.0-1.2 


0.303 


1 

1.8 


6 (6) 


7(7) 




1.2-1.4 


0.305 


1.7 


7(7) 


8 (8) 




1.4-1.6 


0.322 


1.7 


5 (4) 


7(6) 




1.6-1.8 


0.301 


1.7 


6 (5) 


7 (7) 




1.8-2 


0.248 


1.9 


6 (5) 


7 (7) 




2-2.2 


0.159 


2.3 


7(5) 


8 (7) 




L.L—LA 


U.Uc>i4 




(Dj 


( n\ 
9(7) 




0-2.4 


3.08 


0.5 


6(6) 


8(7) 




0-0.4 


0.158 


2.4 


8(7) 


9(9) 




0.4-0.8 


0.170 


2.6 


8(7) 


9(9) 


Y(2S) 


0.8-1.2 


0.179 


2.8 


7(7) 


9(9) 




1.2-1.6 


0.185 


2.7 


6(5) 


8(7) 




1.6-2 


0.157 


2.9 


6(6) 


8(7) 




2-2.4 


0.0577 


7.6 


7(6) 


11 (11) 




0-2.4 


0.907 


1.2 


6(5) 


7(7) 




0-0.4 


0.0858 


3.7 


8(8) 


10(10) 


Y(3S) 


0.4-0.8 


0.0946 


4.0 


8(7) 


10 (9) 




0.8-1.4 


0.141 


3.8 


7(7) 


9(9) 




1.4-2 


0.134 


4.0 


7(7) 


9(9) 




2-2.4 


0.0321 


10.3 


7(5) 


13(12) 




0-2.4 


0.487 


2.0 


6(6) 


8(7) 



Table A. 11: The ratio of Y(3S) /Y(1S) fiducial cross sections, and its px dependence, integrated 
over the rapidity range < 2.4, with statistical and systematic uncertainties combined in 
quadrature. The numbers in parentheses are negative variations. 



px(GeV/c) 



Y(3S)/Y(1S) 



0-2 0.129 ± 0.009 (0.008) ± 0.017 (0.020) 

2-5 0.141 ± 0.008 (0.008) ± 0.018 (0.016) 

5-8 0.154 ± 0.009 (0.009) ± 0.016 (0.016) 

8-10 0.172 ± 0.011 (0.011) ± 0.022 (0.014) 

10-13 0.197 ± 0.012 (0.012) ± 0.022 (0.016) 

13-16 0.251 ± 0.019 (0.018) ± 0.037 (0.021) 

16-18 0.308 ± 0.035 (0.033) ± 0.031 (0.037) 

18-22 0.287 ± 0.029 (0.027) ± 0.023 (0.037) 

22-38 0.345 ± 0.042 (0.038) ± 0.055 (0.051) 

0-38 0.158 ± 0.004 (0.004) ± 0.016 (0.016) 
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Table A. 12: The ratio of Y(2S)/Y(1S) fiducial cross sections, and its pj dependence, integrated 
over the rapidity range \y^\ < 2.4, with statistical and systematic uncertainties combined in 
quadrature. The numbers in parentheses are negative variations. 

PT(GeV/c) Y(2S)/Y(1S) 

0-2 0.258 ± 0.011 (0.011) ± 0.034 (0.032) 

2-5 0.283 ± 0.010 (0.010) ± 0.035 (0.031) 

5-8 0.302 ± 0.012 (0.012) ± 0.031 (0.028) 

8-10 0.289 ± 0.014 (0.013) ± 0.030 (0.022) 

10-13 0.339 ± 0.016 (0.016) ± 0.032 (0.029) 

13-16 0.376 ± 0.023 (0.022) ± 0.053 (0.034) 

16-18 0.481 ± 0.045 (0.042) ± 0.046 (0.053) 

18-22 0.454 ± 0.037 (0.035) ± 0.040 (0.053) 

22-38 0.511 ± 0.051 (0.047) ± 0.083 (0.046) 

0-38 0.296 ± 0.005 (0.005) ± 0.031 (0.029) 



Table A.13: The ratio of Y(3S)/Y(2S) fiducial cross sections, and its pj dependence, integrated 
over the rapidity range |y^| < 2.4, with statistical and systematic uncertainties combined in 
quadrature. The numbers in parentheses are negative variations. 



pj{GeV/c) Y(3S)/Y(2S) 



0- 


-2 


0.500 ± 0.043 (0.040) ± 0.062 (0.072) 


2- 


-5 


0.497 ± 0.036 (0.034) ± 0.063 (0.057) 


5- 


-8 


0.510 ± 0.039 (0.037) ± 0.053 (0.051) 


8- 


10 


0.595 ± 0.049 (0.046) ± 0.079 (0.062) 


10- 


-13 


0.581 ± 0.046 (0.043) ± 0.067 (0.049) 


13- 


-16 


0.669 ± 0.061 (0.057) ± 0.066 (0.054) 


16- 


-18 


0.639 ± 0.089 (0.079) ± 0.083 (0.077) 


18- 


-22 


0.631 ± 0.077 (0.070) ± 0.057 (0.065) 


22- 


-38 


0.675 ± 0.096 (0.085) ± 0.068 (0.101) 


0- 


38 


0.534 ± 0.017 (0.016) ± 0.054 (0.051) 
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Table A.14: The product of the Y(1S) acceptance-corrected production cross sections, a, inte- 
grated and differential in -p^, and the dimuon branching fraction, B, measured for four polar- 
ization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for = 1 and 
\q = —1, integrated over the rapidity range \y^\ < 2.4. The statistical uncertainty (stat.), the 
sum of the systematic uncertainties in quadrature (Egyst.), and the total uncertainty {Aa; includ- 
ing stat., Esyst.' the imcertainty in the integrated luminosity) are in percent. For the four 
polarization scenarios the fractional change to the central value of the cross section relative to 
the unpolarized value is given in percent. The numbers in parentheses are negative variations. 



ly^l < 2.4 
Y(1S) 





rr . K (rth^ 
V Lj yiiuj 


stat. 

a- 


Iwsyst. 


cr 




HX 


1 — —1 
) — J- 


A/, — 1 
Ag — i 


CS 


— J- 


0-0.5 


0.0859 


5.4 


8 (7) 


11 fio) 


+19 




—24 


+21 




—26 


0.5-1 


0.263 


3.3 


8(7) 


9(9) 


+22 




-24 


+23 




-25 


1-1.5 


0.374 


2.6 


8(8) 


9(9) 


+23 




-22 


+26 




-24 


1.5-2 


0.505 


2.4 


9(8) 


10(9) 


+16 




-27 


+19 




-29 


2-3 


1.16 


1.6 


8(10) 


9(11) 


+20 




-25 


+26 




-27 


3-4 


1.21 


2.1 


7(6) 


9(8) 


+20 




-26 


+21 




-26 


4-5 


1.084 


2.1 


7(6) 


8(8) 


+21 




-25 


+21 




-24 


5-6 


0.879 


1.9 


7(9) 


8(10) 


+19 




-25 


+16 




-23 


6-7 


0.680 


2.6 


6(6) 


8(7) 


+18 




-24 


+13 




-23 


7-8 


0.556 


2.0 


6(5) 


7(7) 


+20 




-22 


+12 




-13 


8-9 


0.419 


2.2 


5(5) 


7(7) 


+19 




-22 


+9 




-9 


9-10 


0.331 


2.3 


5(4) 


7(6) 


+19 




-21 


+7 




-4 


10-11 


0.238 


2.5 


5(4) 


7(6) 


+17 




-20 


+4 




-1 


11-12 


0.179 


2.9 


5(4) 


7(6) 


+18 




-19 


+6 




+0 


12-13 


0.145 


2.9 


5(4) 


7(7) 


+16 




-20 


+2 




+2 


13-14 


0.0990 


3.2 


4(5) 


7(7) 


+14 




-21 


-0 




+3 


14-15 


0.0750 


3.6 


5(5) 


8(7) 


+15 




-20 


-1 




+5 


15-16 


0.0595 


3.8 


5(5) 


7(7) 


+13 




-19 


-1 




+5 


16-18 


0.0732 


3.4 


5(5) 


7(7) 


+12 




-20 


-3 




+6 


18-20 


0.0500 


3.8 


5(4) 


7(7) 


+13 




-18 


-2 




+7 


20-22 


0.0302 


5.1 


5(4) 


8(8) 


+11 




-18 


-3 




+6 


22-25 


0.0237 


5.6 


5(4) 


8(8) 


+9 




-18 


-3 




+4 


25-30 


0.0205 


6.0 


5(4) 


9(8) 


+10 




-16 


-1 




+4 


30-50 


0.0123 


7.4 


6(6) 


10(10) 


+5 




-17 


-7 




+3 


0-50 


8.55 


0.6 


7(6) 


8(7) 


+19 




-24 


+16 




-19 
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Table A.15: The product of the Y(2S) acceptance-corrected production cross sections, a, inte- 
grated and differential in p^, and the dimuon branching fraction, B, measured for four polar- 
ization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for = 1 and 
Ag = —1, integrated over the rapidity range \y^\ < 2.4. The statistical uncertainty (stat.), the 
sum of the systematic uncertainties in quadrature (Egyst.), and the total uncertainty (Atr; includ- 
ing stat., Esyst.' the imcertainty in the integrated luminosity) are in percent. For the four 
polarization scenarios the fractional change to the central value of the cross section relative to 
the unpolarized value is given in percent. The numbers in parentheses are negative variations. 



ly^l <2.4 

Y(2S) 



pT(GeV/c) 


(T-B (nb) 


stat. 

cr 


I-syst. 




A0 = 1 


HX 


) = -l 


A0 = 1 


CS 

At 


) = -1 


0-1 


0.0829 


5.9 


9(8) 


11 (11) 


+18 




-23 


+19 




-24 


1-2.5 


0.331 


3.3 


11 (10) 


12(11) 


+14 




-23 


+18 




-25 


2.5^ 


0.409 


3.1 


9(8) 


10 (9) 


+22 




-24 


+21 




-26 


4-5.5 


0.362 


3.3 


8(7) 


9(9) 


+18 




-22 


+16 




-24 


5.5-7 


0.286 


3.6 


7(6) 


9(8) 


+17 




-23 


+15 




-20 


7-8.5 


0.212 


3.9 


7(7) 


9(9) 


+21 




-20 


+15 




-12 


8.5-10 


0.146 


4.0 


6(6) 


9(8) 


+20 




-18 


+11 




-5 


10-11.5 


0.1123 


4.1 


6(6) 


9(8) 


+19 




-18 


+8 




-0 


11.5-13 


0.0765 


4.6 


5(5) 


8(8) 


+17 




-18 


+6 




+1 


13-14.5 


0.0519 


5.1 


5(5) 


8(8) 


+14 




-20 


+0 




+3 


14.5-16 


0.0376 


5.7 


5(7) 


9(10) 


+14 




-19 


+1 




+3 


16-18 


0.0373 


5.3 


6(5) 


9(8) 


+13 




-18 


+1 




+3 


18-19.5 


0.0159 


7.4 


5(4) 


10 (9) 


+13 




-18 


-3 




+9 


19.5-22 


0.0204 


6.6 


5(5) 


9(9) 


+11 




-16 


-1 




+5 


22-26 


0.0158 


7.2 


5(4) 


10 (9) 


+11 




-16 


-3 




+7 


26^2 


0.0126 


7.7 


6(5) 


10(10) 


+10 




-15 


-4 




+9 


0^2 


2.21 


1.2 


7(6) 


8(7) 


+14 




-24 


+13 




-20 
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Table A. 16: The product of the Y(3S) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the dimuon branching fraction, B, measured for four polar- 
ization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for = 1 and 
Xq = —1, integrated over the rapidity range \y^\ < lA. The statistical uncertainty (stat.), the 
sum of the systematic uncertainties in quadrature (Egyst.), and the total uncertainty (Atr; includ- 
ing stat., I^syst./ and the uncertainty in the integrated luminosity) are in percent. For the four 
polarization scenarios the fractional change to the central value of the cross section relative to 
the unpolarized value is given in percent. The numbers in parentheses are negative variations. 



ly^l < 2.4 
Y(3S) 



pT(GeV/c) 


a ■ B (nb) 


stat. 

a 




Aa 


HX 

\e = l h 


) = -1 


A0 = 1 


CS 

Ke = -1 


0-2.5 


0.203 


5.3 


8(8) 


11 (10) 


+17 


-21 


+20 


-25 


2.5-5 


0.287 


4.5 


10 (11) 


12 (12) 


+17 


-22 


+20 


-25 


5-7.5 


0.227 


4.6 


9(8) 


11 (10) 


+16 


-22 


+20 


-22 


7.5-10 


0.157 


4.8 


11 (10) 


12(12) 


+23 


-16 


+16 


-5 


10-13 


0.113 


4.3 


7(5) 


9(8) 


+20 


-15 


+12 


-1 


13-16 


0.0617 


5.0 


5(5) 


8(8) 


+14 


-17 


+4 


+1 


16-18 


0.0227 


7.4 


6(5) 


10(10) 


+12 


-18 


-1 


+4 


18-22 


0.0229 


7.0 


7(6) 


10(10) 


+12 


-17 


-1 


+6 


22-38 


0.0185 


7.6 


13(13) 


15(15) 


+10 


-15 


-0 


+6 


0-38 


1.11 


2.0 


9(8) 


10(9) 


+16 


-21 


+14 


-17 
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Table A.17: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
= 1 and Ag = —1, integrated over the rapidity range < 0.4. The statistical uncertainty 
(stat.), the sum of the systematic imcertainties in quadrature (Egyst.)/ ^rid the total uncertainty 
(Act; including stat., I^gygt./ arid the uncertainty in the integrated luminosity) are in percent. For 
the four polarization scenarios the fractional change to the central value of the cross section 
relative to the unpolarized value is given in percent. The niuxibers in parentheses are negative 
variations. 



ly^l < 0-4 















HX 






CS 




PT(GeV/c) 


a ■ B (nb) 


stat. 

a 




Aa 


Ae = l At 


) = -1 


Ae = l 


Ae = -1 




0-2 


0.216 


3.4 


14(13) 


15 (14) 


+10 


-16 


+21 


-27 




2-4 


0.387 


3.2 


11 (10) 


12 (12) 


+15 


-22 


+20 


-26 




4r-6 


0.355 


3.0 


10(10) 


11 (11) 


+21 


-27 


+16 


-22 


Y(1S) 


6-8 


0.224 


3.2 


9(9) 


10 (10) 


+17 


-30 


+3 


-17 


8-11 


0.190 


2.9 


7(7) 


8(8) 


+22 


-27 


-0 


+1 




11-15 


0.0914 


3.3 


5(5) 


7(7) 


+20 


-25 


-4 


+10 




15-50 


0.0509 


3.7 


4(4) 


7(7) 


+17 


-22 


-5 


+13 




0-3 


0.084 


6.6 


16(15) 


18 (17) 


+22 


-13 


+20 


-25 




3-7 


0.143 


4.9 


13(16) 


15 (17) 


+27 


-16 


+15 


-19 


Y(2S) 


7-11 


0.0813 


5.3 


9(8) 


11 (11) 


+21 


-26 


+2 


-3 


11-15 


0.0323 


6.1 


6(5) 


9(9) 


+20 


-25 


-4 


+9 




15^2 


0.0243 


5.7 


5(4) 


8(8) 


+16 


-21 


-5 


+12 




0-7 


0.100 


6.5 


15 (14) 


17(16) 


+12 


-16 


+16 


-22 


Y(3S) 


7-12 


0.0487 


7.5 


11 (11) 


14 (14) 


+20 


-25 


+2 


-3 


12-38 


0.0291 


5.9 


5(5) 


9(9) 


+17 


-23 


-5 


+11 
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Table A.18: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
\g = 1 and Ag = —1, integrated over the rapidity range 0.4 < < 0.8. The statistical 
uncertainty (stat.), the sum of the systematic uncertainties in quadrature (Egyst.)' the total 
uncertainty (Au; including stat., I^gyst./ arid the uncertainty in the integrated luminosity) are in 
percent. For the four polarization scenarios the fractional change to the central value of the 
cross section relative to the unpolarized value is given in percent. The numbers in parentheses 
are negative variations. 



0.4 < ly^l < 0.8 















HX 






CS 




PT(GeV/c) 


a ■ B (nb) 


stat. 

(7 




ha 


Ae = 1 Af 


) = -1 


\e = l 


Ae = -1 




0-2 


0.220 


4.5 


13 (12) 


14(13) 


+20 


-25 


+21 


-26 




2-4 


0.409 


3.0 


10 (9) 


11 (10) 


+21 


-25 


+22 


-25 




4r-6 


0.367 


2.9 


8(8) 


9(9) 


+20 


-25 


+16 


-21 


Y(1S) 


6-8 


0.231 


3.3 


7(6) 


9(8) 


+13 


-27 


+2 


-15 


8-11 


0.180 


2.9 


6(5) 


8(7) 


+19 


-23 


-0 


+2 




11-15 


0.0915 


3.3 


6(6) 


8(8) 


+18 


-23 


-4 


+10 




15-50 


0.0492 


3.8 


4(4) 


7(7) 


+15 


-20 


-5 


+13 




0-3 


0.088 


7.3 


14(13) 


16 (15) 


+18 


-23 


+19 


-24 




3-7 


0.172 


4.5 


10 (9) 


11 (11) 


+17 


-23 


+15 


-20 


Y(2S) 


7-11 


0.0811 


5.2 


6(6) 


9(9) 


+17 


-23 


+1 


-2 


11-15 


0.0293 


6.6 


5(5) 


9(9) 


+17 


-22 


-4 


+9 




15^2 


0.0227 


5.9 


4(4) 


8(8) 


+15 


-20 


-5 


+12 




0-7 


0.115 


6.7 


14 (13) 


16 (15) 


+16 


-22 


+16 


-21 


Y(3S) 


7-12 


0.0558 


6.8 


7(6) 


11 (10) 


+16 


-22 


+1 


-2 


12-38 


0.0257 


6.5 


5(5) 


9(9) 


+15 


-21 


-4 


+10 
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Table A.19: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
\g = 1 and \g = —1, integrated over the rapidity range 0.8 < |y^| < 1.2 The statistical un- 
certainty (stat.), the sum of the systematic uncertainties in quadrature (Zgyst.)' the total 
uncertainty (Au; including stat., I^gyst./ arid the uncertainty in the integrated luminosity) are in 
percent. For the four polarization scenarios the fractional change to the central value of the 
cross section relative to the unpolarized value is given in percent. The numbers in parentheses 
are negative variations. 



0.8 < ly^l < 1.2 















HX 






CS 




PT(GeV/c) 


a ■ B (nb) 


stat. 

(7 




ha 


Ae = 1 At 


) = -1 


\e = l 


Ae = -1 




0-2 


0.198 


3.7 


11 (9) 


12 (10) 


+20 


-25 


+20 


-24 




2-4 


0.426 


2.7 


9(9) 


10 (10) 


+20 


-26 


+19 


-25 




4r-6 


0.331 


2.9 


7(6) 


8(8) 


+17 


-25 


+13 


-21 


Y(1S) 


6-8 


0.217 


3.4 


7(7) 


9(8) 


+22 


-18 


+12 


-7 


8-11 


0.174 


2.9 


5(5) 


7(7) 


+15 


-21 


-1 


+2 




11-15 


0.0879 


3.4 


6(5) 


8(8) 


+15 


-21 


-4 


+9 




15-50 


0.0482 


3.8 


5(5) 


8(8) 


+13 


-19 


-5 


+12 




0-3 


0.098 


6.2 


14(13) 


15 (15) 


+18 


-23 


+17 


-23 




3-7 


0.155 


5.3 


11 (8) 


13 (10) 


+19 


-22 


+16 


-19 


Y(2S) 


7-11 


0.0803 


5.3 


8(7) 


10 (10) 


+14 


-20 


+1 


-2 


11-15 


0.0309 


6.2 


6(6) 


10 (9) 


+14 


-20 


-4 


+9 




15^2 


0.0214 


6.2 


6(6) 


10 (9) 


+13 


-18 


-5 


+11 




0-7 


0.109 


7.7 


15 (14) 


17(17) 


+16 


-22 


+14 


-20 


Y(3S) 


7-12 


0.0579 


6.9 


9(9) 


12 (12) 


+13 


-19 


+2 


-3 


12-38 


0.0250 


6.8 


6(6) 


10 (10) 


+13 


-19 


-4 


+10 
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Table A.20: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
Xq = 1 and = —1, integrated over the rapidity range 1.2 < |y^| < 1.6 The statistical un- 
certainty (stat.), the sum of the systematic uncertainties in quadrature (Egyst.)/ aJtid the total 
uncertainty (Au; including stat., I^syst./ ^rid the uncertainty in the integrated luminosity) are in 
percent. For the four polarization scenarios the fractional change to the central value of the 
cross section relative to the unpolarized value is given in percent. The niunbers in parentheses 
are negative variations. 



1.2 < ly^l < 1.6 

















HX 




CS 




pT(GeV/c) 


a ■ B (nb) 


stat. 

(7 


X^syst. 


ha 


Ae = 1 


Ae = -1 


Xe = l 


Xe = -1 




0-2 


0.203 


3.5 


10 (8) 


11 (9) 


+19 


-23 


+21 


-22 




2-4 


0.416 


2.4 


7(6) 


8(7) 


+20 


-26 


+18 


-25 




^6 


0.315 


2.8 


6(5) 


8(7) 


+20 


-24 


+17 


-21 


Y(1S) 


6-8 


0.215 


3.0 


5(6) 


7(7) 


+17 


-22 


+11 


-15 


8-11 


0.165 


3.0 


5(5) 


7(7) 


+13 


-19 


+2 


-4 




11-15 


0.0868 


3.3 


5(5) 


7(7) 


+13 


-19 


-3 


+6 




15-50 


0.0485 


3.8 


4(4) 


7(7) 


+12 


-18 


-5 


+11 




0-3 


0.099 


6.2 


11 (11) 


13 (13) 


+18 


-24 


+18 


-23 




3-7 


0.161 


4.5 


7(6) 


9(9) 


+18 


-23 


+17 


-22 


Y(2S) 


7-11 


0.0758 


5.5 


7(6) 


10(9) 


+14 


-19 


+6 


-9 


11-15 


0.0305 


6.5 


5(5) 


9(9) 


+12 


-18 


-2 


+5 




15-42 


0.0219 


6.4 


5(4) 


9(9) 


+11 


-17 


-4 


+10 




0-7 


0.113 


6.8 


7(7) 


11 (10) 


+17 


-23 


+16 


-22 


Y(3S) 


7-12 


0.0480 


7.7 


10(10) 


13 (13) 


+13 


-18 


+5 


-8 


12-38 


0.0276 


6.7 


8(7) 


11 (11) 


+11 


-17 


-3 


+7 
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Table A.21: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
Xq = 1 and = —1, integrated over the rapidity range 1.6 < |y^| < 2.0. The statistical 
imcertainty (stat.), the sum of the systematic uncertainties in quadrature (Egyst.)/ arid the total 
uncertainty (Act; including stat., I^syst.' ^rid the uncertainty in the integrated luminosity) are in 
percent. For the four polarization scenarios the fractional change to the central value of the 
cross section relative to the unpolarized value is given in percent. The niunbers in parentheses 
are negative variations. 



1.6 < ly^l < 2.0 















HX 




CS 






pT(GeV/c) 


a ■ B (nb) 


stat. 

a 


X^syst. 

a 


(7 


A(5 = 1 Af 


1 = -1 


\e = l h 


= -1 




0-2 


0.211 


3.1 


7(6) 


9(8) 


+22 


-26 


+23 


-27 




2-4 


0.401 


2.5 


6(4) 


7(6) 


+22 


-26 


+24 


-27 




^6 


0.330 


2.6 


5(4) 


7(6) 


+20 


-25 


+23 


-27 


Y(1S) 


6-8 


0.191 


3.2 


4(4) 


7(6) 


+18 


-22 


+20 


-24 


8-11 


0.151 


3.4 


5(4) 


7(7) 


+15 


-20 


+15 


-19 




11-15 


0.0746 


4.1 


5(5) 


8(7) 


+10 


-16 


+5 


-9 




15-50 


0.0465 


4.4 


6(6) 


9(8) 


+9 


-14 


-1 


+3 




0-3 


0.091 


6.8 


8(8) 


11 (11) 


+21 


-25 


+22 


-27 




3-7 


0.175 


4.4 


7(6) 


9(9) 


+19 


-24 


+23 


-27 


Y(2S) 


7-11 


0.0608 


7.0 


7(6) 


10 (10) 


+14 


-20 


+17 


-22 


11-15 


0.0263 


8.3 


5(5) 


11 (10) 


+10 


-15 


+7 


-11 




15-42 


0.0205 


8.1 


6(5) 


11 (11) 


+8 


-13 


-0 


+1 




0-7 


0.117 


7.2 


9(9) 


13 (12) 


+19 


-23 


+23 


-26 


Y(3S) 


7-12 


0.0341 


10.6 


18 (18) 


22(21) 


+13 


-19 


+16 


-21 


12-38 


0.0243 


8.8 


6(6) 


11 (11) 


+8 


-13 


+4 


-6 
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Table A.22: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
Xq = 1 and \e = —1, integrated over the rapidity range 2.0 < |y^| < 2.4 . The statistical 
imcertainty (stat.), the sum of the systematic uncertainties in quadrature (Egyst.)/ arid the total 
uncertainty (Act; including stat., I^syst.' ^rid the uncertainty in the integrated luminosity) are in 
percent. For the four polarization scenarios the fractional change to the central value of the 
cross section relative to the unpolarized value is given in percent. The niunbers in parentheses 
are negative variations. 



2.0 < ly^l < 2.4 















HX 




CS 






pT(GeV/c) 


a ■ B (nb) 


stat. 

a 


X^syst. 


ha 


A0 = 1 h 


) = -1 




) = -1 




0-2 


0.135 


5.3 


8(7) 


11(9) 


+30 


-31 


+32 


-31 




2-4 


0.338 


2.8 


7(6) 


9(8) 


+23 


-28 


+30 


-32 




^6 


0.289 


4.1 


8(6) 


10(9) 


+27 


-21 


+38 


-28 


Y(1S) 


6-8 


0.150 


5.8 


7(6) 


10(9) 


+23 


-21 


+34 


-29 


8-11 


0.126 


5.8 


7(5) 


10(9) 


+18 


-20 


+31 


-29 




11-15 


0.0667 


7.2 


8(7) 


12(11) 


+14 


-18 


+28 


-29 




15-50 


0.0356 


9.1 


9(9) 


14(13) 


+6 


-13 


+18 


-24 




0-3 


0.0475 


13.2 


9(8) 


16 (16) 


+27 


-30 


+30 


-32 




3-7 


0.148 


7.5 


7(5) 


11 (10) 


+20 


-25 


+31 


-32 


Y(2S) 


7-11 


0.0636 


10.9 


8(7) 


14 (13) 


+16 


-22 


+29 


-31 


11-15 


0.0224 


14.7 


9(8) 


18 (17) 


+12 


-18 


+26 


-29 




15-42 


0.0137 


20.4 


11 (11) 


24 (23) 


+9 


-14 


+20 


-24 




0-7 


0.105 


11.2 


8(6) 


14(13) 


+22 


-28 


+30 


-32 


Y(3S) 


7-12 


0.0366 


19.3 


11 (9) 


22(22) 


+14 


-20 


+29 


-31 


12-38 


0.0231 


15.8 


9(8) 


19 (18) 


+9 


-14 


+24 


-27 
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Table A.23: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, measured for 
four polarization scenarios, in the helicity frame (HX) and Colins-Soper (CS) frame, each for 
= 1 and = —1, integrated over the px range pj < 50GeV/c. The statistical uncertainty 
(stat.), the sum of the systematic imcertainties in quadrature (Egyst.), and the total uncertainty 
(Au; including stat., I]syst./ arid the uncertainty in the integrated luminosity) are in percent. For 
the four polarization scenarios the fractional change to the central value of the cross section 
relative to the unpolarized value is given in percent. The niunbers in parentheses are negative 
variations. 



p| <50GeV/c 

















rlA 








\y\ 


a ■ B (nb) 


stat. 

a 


(7 


a 


A0 = 1 


A0 = -1 


Af, = l 


A0 = -1 




0-0.2 


0.770 


1.7 


10 (9) 


11 (10) 


+27 


-16 


+22 


-9 




0.2-0.4 


0.777 


1.7 


9(9) 


10(10) 


+19 


-25 


+12 


-16 




0.4-0.6 


0.770 


1.7 


9(9) 


10 (10) 


+20 


-18 


+13 


-15 


Y(1S) 


0.6-0.8 


0.795 


1.7 


8(8) 


9(9) 


+20 


-24 


+13 


-15 




0.8-1.0 


0.761 


1.7 


9(8) 


10 (9) 


+18 


-24 


+12 


-15 




1.0-1.2 


0.738 


1.8 


6(5) 


7(7) 


+18 


-24 


+11 


-16 




1. 2-1 A 


0.714 


1.7 


7(7) 


8(8) 


+18 


-23 


+12 


-16 




1.4-1.6 


0.738 


1.7 


5(4) 


7(6) 


+18 


-23 


+14 


-18 




1.6-1.8 


0.721 


1.7 


6(5) 


7(7) 


+19 


-23 


+18 


-21 




1.8-2 


0.708 


1.9 


6(5) 


7(6) 


+20 


-24 


+23 


-26 




2-2.2 


0.636 


2.3 


6(5) 


8(7) 


+20 


-24 


+28 


-30 




2.2-2.4 


0.528 


4.2 


7(5) 


9(8) 


+22 


-26 


+32 


-32 




0-2.4 


8.66 


0.6 


7(6) 


8(7) 


+19 


-24 


+16 


-19 




0-0.4 


0.376 


2.5 


10(9) 


11(11) 


+16 


-21 


+10 


-13 




0.4-0.8 


0.400 


2.5 


9(8) 


10 (9) 


+17 


-23 


+10 


-13 


Y(2S) 


0.8-1.2 


0.380 


2.8 


9(8) 


10 (9) 


+16 


-22 


+10 


-13 




1.2-1.6 


0.382 


2.7 


7(6) 


8(8) 


+16 


-22 


+12 


-16 




1.6-2 


0.380 


2.9 


7(6) 


8(8) 


+17 


-22 


+19 


-23 




2-2.4 


0.295 


5.1 


7(5) 


9(8) 


+19 


-24 


+29 


-31 




0-2.4 


2.21 


1.3 


8(7) 


9(8) 


+26 


-16 


+23 


-12 




0-0.4 


0.185 


3.9 


11 (10) 


12 (12) 


+15 


-20 


+9 


-11 


Y(3S) 


0.4-0.8 


0.202 


3.9 


10(9) 


11(11) 


+16 


-22 


+9 


-11 




0.8-1.4 


0.282 


3.8 


9(9) 


11 (10) 


+15 


-21 


+8 


-11 




1.4-2 


0.290 


4.0 


7(7) 


9(9) 


+16 


-21 


+16 


-20 




2-2.4 


0.164 


8.3 


7(6) 


12(11) 


+18 


-23 


+29 


-31 




0-2.4 


1.12 


2.1 


8(7) 


10(8) 


+16 


-21 


+14 
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Table A.24: The ratio of Y(3S)/Y(1S) acceptance-corrected cross sections, and its pj depen- 
dence, integrated over the rapidity range < 2.4, for the unpolarized scenario, with statisti- 
cal and systematic uncertainties. The niuxibers in parentheses are negative variations. 



PT(GeV/c) Y(3S)/Y(1S) 

0-2 0.108 ± 0.007 (0.007) ± 0.015 (0.015) 

2-5 0.101 ± 0.006 (0.006) ± 0.015 (0.012) 

5-8 0.122 ± 0.007 (0.007) ± 0.016 (0.014) 

8-10 0.154 ± 0.010 (0.009) ± 0.022 (0.017) 

10-13 0.197 ± 0.012 (0.011) ± 0.020 (0.015) 

13-16 0.258 ± 0.019 (0.018) ± 0.021 (0.019) 

16-18 0.305 ± 0.034 (0.032) ± 0.027 (0.026) 

18-22 0.283 ± 0.030 (0.028) ± 0.022 (0.027) 

22-38 0.323 ± 0.039 (0.036) ± 0.050 (0.089) 

0-38 0.129 ± 0.003 (0.003) ± 0.018 (0.015) 



Table A.25: The ratio of Y(2S)/Y(1S) acceptance-corrected cross sections, and its pj depen- 
dence, integrated over the rapidity range | < 2.4, for the unpolarized scenario, with statisti- 
cal and systematic uncertainties. The numbers in parentheses are negative variations. 



PT(GeV/c) Y(2S)/Y(1S) 

0-2 0.221 ± 0.010 (0.009) ± 0.030 (0.024) 

2-5 0.230 ± 0.008 (0.008) ± 0.028 (0.030) 

5-8 0.266 ± 0.011 (0.010) ± 0.033 (0.025) 

8-10 0.270 ± 0.012 (0.012) ± 0.034 (0.029) 

10-13 0.330 ± 0.016 (0.015) ± 0.029 (0.028) 

13-16 0.372 ± 0.023 (0.022) ± 0.025 (0.028) 

16-18 0.502 ± 0.046 (0.043) ± 0.040 (0.041) 

18-22 0.449 ± 0.039 (0.036) ± 0.032 (0.048) 

22-38 0.466 ± 0.047 (0.044) ± 0.074 (0.121) 

0-38 0.264 ± 0.005 (0.005) ± 0.032 (0.028) 



Table A.26: The ratio of Y(3S)/Y(2S) acceptance-corrected cross sections, and its pj depen- 
dence, integrated over the rapidity range < 2.4, for the unpolarized scenario, with statisti- 
cal and systematic uncertainties. The nimibers in parentheses are negative variations. 



PT(GeV/c) Y(3S)/Y(2S) 



0- 


-2 


0.490 ± 0.042 (0.040) ± 0.067 (0.070) 


2-5 


0.440 ± 0.032 (0.030) ± 0.082 (0.057) 


5- 


-8 


0.459 ± 0.035 (0.033) ± 0.054 (0.053) 


8- 


10 


0.569 ± 0.047 (0.044) ± 0.078 (0.063) 


10- 


-13 


0.597 ± 0.046 (0.043) ± 0.062 (0.044) 


13- 


-16 


0.693 ± 0.063 (0.059) ± 0.060 (0.047) 


16- 


-18 


0.607 ± 0.082 (0.074) ± 0.058 (0.051) 


18- 


-22 


0.630 ± 0.079 (0.071) ± 0.060 (0.044) 


22- 


-38 


0.69 ± 0.10 (0.09) ± 0.15 (0.14) 


0- 


■38 


0.490 ± 0.015 (0.015) ± 0.069 (0.058) 



41 



Table A.27: The product of the Y(nS) acceptance-corrected production cross sections, a, inte- 
grated and differential in pj, and the respective dimuon branching fraction, B, integrated over 
the rapidity range |y^| < 1.2 and the pj range from 10 to 50 GeV/c. The cross sections assume 
the Y(nS) are unpolarized. The statistical uncertainty (stat.), the sum in quadrature of the sys- 
tematic uncertainties excluding the contribution from the polarization uncertainty (Esyst.)/ the 
systematic uncertainties from the polarization (syst.(pol)), and the total uncertainty {Aa; in- 
cluding stat., Esyst.' syst.(pol), and the imcertainty in the integrated luminosity) are in percent. 
The numbers in parentheses are negative variations. 
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